Majorana Fermions in Nodal Superconductors: Gapless Topological Phase 
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We demonstrate that Majorana fermions exist in edges of systems and in a vortex core even 
for superconductors with nodal excitations such as the d-wave pairing state under a particular but 
realistic condition in the case with an anti-symmetric spin-orbit interaction and a nonzero magnetic 
field below the upper critical field. We clarify that the Majorana fermion state is topologically 
protected in spite of the presence of bulk gapless nodal excitations, because of the existence of a 
nontrivial topological number. Our finding drastically enlarges target systems where we can explore 
the Majorana fermion state. 

PACS numbers: 



Introduction - - Majorana fermions realized as vor- 
tex core bound states of superconducting condensates 
have been attracting considerable interest in connec- 
tion with the application to quantum computationjli- 
.3]. Such vortices obey the non-Abelian statistics [4J-19J, 
and because of this distinct feature, they are utilized 
as decoherence-free qubits. The realization of Majo- 
rana bound states has been discussed for the quantum 
Hall effect systems [#tZ|, p + ip superconductors M - 11 1 , 
supcrconductor-topological- insulator interfaces |12l Il3| . 
and s-wave Rashba superconductors 14Hl6| | . The origin 
of Majorana fermions acting as non-Abelian anyons is 
intimately related to the existence of the non-Abelian 
topological order, which yields the fractionalization of 
quasiparticles[l7l |l8j. Generally, topological order is 
characterized by a nontrivial topological number asso- 
ciated with the global structure of the Hilbert space, and 
hence, the existence of a nonzero energy gap which sep- 
arates the topological ground state and non-topological 
excited states stabilizes the topological order. 

In this paper, we propose an example of systems real- 
izing Majorana fermions, which is unusual in the above- 
mentioned sense of topological stability, but ubiquitous 
in real materials: Majorana fermion states realized in su- 
perconductors with nodal excitations such as the d-wave 
pairing state. More precisely, Majorana fermions exist in 
edges of systems and in a vortex core for nodal supercon- 
ductors, provided that there are both an anti-symmetric 
spin-orbit interaction and a nonzero magnetic field below 
the upper critical field H C 2, and that when the magnetic 
field is switched off, the Fermi level is located close to 
odd numbers of time-reversal invariant fc-points in the 
Brillouin zone at which the superconducting gap van- 
ishes because of symmetry requirement. This proposal 
implies that the non-Abelian topological order coexists 
with gapless excitations in our system. One may wonder 
how the topological stability is ensured in the presence 
of non-topological gapless excitations. In fact, the Chern 
number is not well-defined in our nodal system. Never- 
theless, we clarify that the Majorana fermion state is, to 



some extent, stable against interactions with nodal exci- 
tations and with impurities, because of the existence of 
a topological number which is well-defined even for gap- 
less superconductors. Note that many classes of noncen- 
trosymmetric (NCS) superconductors such as CePtaSi, 
CeRhSi3, CeIrSi3, L^PfaB, are known to possess super- 
conducting gap-nodes 19-23]. In these systems, some of 
time-reversal invariant fc-points reside close to the Fermi 
level 2J] . Our finding indicates that if the total num- 
ber of these fc-points is odd, and the superconducting 
gap vanishes (or, at least, becomes sufficiently small) at 
these points, stable Majorana fermion modes appear un- 
der applied magnetic fields. We expect that such Ma- 
jorana fermion states may be realized in large classes of 
NCS superconductor with gap-nodes. 

Majorana fermions in edges and in a vortex core — To 
be concrete, we consider a two-dimensional d-wave super- 
conductor with the Rashba spin-orbit (SO)interaction, 
though the following argument is basically applicable to 
any NCS nodal superconductors. The Hamiltonian is 
given by H = \ J2 k V'fe'^WV'fe, with 
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where ip' k = (c' h 

cosky) — [i, g k = 2A(sinfc y , — smk x , 0), fc = (k x , k y ), and 
a = (a x , o-y, a z ) the Pauli matrices, h = [1bH z is a Zee- 
man magnetic field. The gap function is the d x 2_ y 2-w&ve 
type, Afc = A()(cosk x — cosk v ), or the d^-wave type, 
Afc = Ao sin k x sin k y . We neglect the orbital effect of 
the magnetic field for a while since it does not change 
our results qualitatively, as long as H z < H c i. 

We, first, demonstrate that there is a gapless chiral 
Majorana fermion mode on the edge of the system. For 
this purpose, we numerically calculate the energy spec- 
trum of the system with the open boundary condition 
imposed for the rc-axis, and the periodic boundary con- 
dition for the the y-axis. The results are shown in FIGH] 
In the case of the d x 2_ y 2-w&ve pairing, when the condi- 
tion — At — /i < h < — /i is satisfied, a gapless edge mode 



appears for k y ~ 0. (FIGQJA)) Note that this edge mode for the d xy -wa.ve pairing, because of the topological ar- 



is isolated from continuum of gapless excitations from the 
gap nodes with the finite Fermi momentum. Because of 
the particle-hole symmetry of the Hamiltonian, the ex- 
istence of one zero energy mode implies that the zero 
energy edge mode for k y ~ is a Majorana fermion. 
The above condition implies that when the Fermi level 
crosses fc-points close to time-reversal invariant points in 
the Brillouin zone in the absence of the magnetic field, 
say, the Y point; i.e. fi = —At, the Majorana edge mode 
appears for any nonzero H z below H c i- This property 
makes a sharp contrast to the s-wave pairing state con- 
sidered in refs. l44l6|. for which a large magnetic field 
satisfying h > A is required to realize the topological or- 
der and Majorana fermions. This is because that the gap 
function of the d-wave pairing vanishes at fc ~ 0, fulfilling 
the condition h > Ak^o ~ 0. As a result, the realization 
of Majorana fermions in the d-wave pairing state is much 
more feasible than that in the s-wave pairing state, which 
may be seriously affected by the orbital depairing effect 
due to the large magnetic field h > A. 
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FIG. 1: Energy spectra for systems with open boundaries for 
the i-direction and the periodic boundary condition for the 
y-direction. fj, - -At, X = 0.5£, A = t and h = 2t. (A) (B) 
d x 2_ y 2 wave pairing. (D) (E) d xy wave pairing. The distance 
between two edges is L = 90 [(A) (D)] and L = 30 [(B) (E)]. 
In (A) and (D), chiral gapless edge modes at x = and x = L 
are depicted, respectively, in green and red curves. (C) (F) 
Chiral edge modes counter-propagating on two opposite edges 
for d x 2_ y 2 wave pairing (C) and d xy wave pairing (F). 



In the case of the d^-wave pairing, the zero energy 
edge mode at k y ~ merges into nodal excitations, 
and thus, it is difficult to identify the Majorana mode 
from the numerical result (FIGQJD)). However, we de- 
duce that there is still a Majorana mode at k y ~ also 



gument which will be presented later, as long as L is 
sufficiently large, (also see the argument on the stability 
below.) 

The existence of the chiral Majorana fermion edge 
mode implies that there is a Majorana fermion mode in 
a vortex core of the superconducting condensate when 
the vorticity n is odd. In the case of the d-wave pairing, 
the analysis of the vortex core state is cumbersome, in 
contrast to the s-wave pairing or the p + ip-wave pair- 
ing states, since the gap function of the d-wave state 
is not an eigen state of the orbital angular momentum, 
and, moreover, there is no truly localized bound state 
in a vortex core because of interactions with delocalized 
nodal excitations [25J, [26|. However as in the case of the 
s-wave Rashba superconductor FJ, [27[ , we can construct 
the Majorana zero energy mode from quasiparticles with 
k ~ at least in some parameter regions. The eigen 
function for this zero energy state is (f> T = (u-f, u±,ul,uV) 

with u-f = ie l ~ ~ e f(r), uj. = — ie I_ 2~ f(r) and f(r) = 
i/-^-e _ 2X r for large r. Here n is odd. In the d-wave 
pairing state, in addition to this zero energy mode, there 
are four gapless extended states outside of the vortex core 
which stem from the four gap-nodes[28ll29|. Since the to- 
tal number of the zero energy mode is odd, one Majorana 
mode survives. Thus, we have a zero energy Majorana 
fermion mode in the vortex core. 

Stability of Majorana fermions — The next important 
question is whether Majorana fermions found above are 
stable or not against weak perturbations such as im- 
purities even in the presence of gapless nodal excita- 
tions. As will be shown later, there is indeed a topo- 
logical protection mechanism in spite of the existence of 
bulk gapless excitations. Before discussing the topologi- 
cal mechanism, however, we here present a heuristic ar- 
gument on this issue to grasp an intuitive physical pic- 
ture. Generally, impurity scattering affects the super- 
conducting state with gap-nodes. We consider only the 
case that the density of impurities is sufficiently small 
so that the superconducting gap is not much reduced. 
We, first, consider the case of the d x 2 _ y 2-wave pairing. 
In a semi-infinite system with an open boundary, there 
is only one chiral Majorana edge mode, while there are 
four gapless modes which stem from four nodes of the 
d 2 .2_ l ,2-wave superconducting gap. To generate an en- 
ergy gap in the Majorana spectrum, we need even num- 
bers of Majorana modes which are paired into complex 
fermions. Thus, for this geometry, the chiral Majorana 
edge mode is stable against interactions with nodal exci- 
tations, and also against impurity scattering. However, 
this argument is not applicable to the case with two open 
boundaries at the opposite sides of the system. In this 
case, two counter-propagating chiral Majorana modes re- 
side in the two opposite edges, as depicted in FIGsfljC). 
Interactions between bulk gapless nodal excitations and 



two chiral Majorana modes may give rise to long-range 
tunneling between two Majorana modes. We note that 
such long-range tunneling via nodal excitations does not 
occur in a clean system, because of the mismatch of the 
Fermi momenta of nodal excitations Uf 7^ and that of 
the chiral Majorana modes with which kp ~ for the 
d a .2_ y 2-wave pairing. When there are impurity poten- 
tials, a Majorana mode and nodal excitations on an edge 
can be hybridized via impurity scattering, leading to the 
long-range tunneling of the Majorana fcrmions in two op- 
posite edges: W tun = t(r - r' o )ij(0,y o )j(L,y' o ), where 
^(x,y) is a Majorana fermion operator, rg = (0,2/o) 
and r' = (L,y' ) are the positions of impurities and 
the tunneling amplitude t(r — r ) ~ l/\r — r' \ for a 
large |r — r |. We introduce a complex fermion opera- 
tor, a(y) = \ [7(0,?/) + ij{L, -y + y + y' )]. Then, the 
Hamiltonian for two chiral Majorana edge states can be 
rewritten into that of the ID chiral Dirac fermion[30j, 
Hodge = —ivjdyar(y)d y a(y). The long-range tunneling 
term is also expressed as Htun — £0 (2a 1, (2/0)0(2/0) — !)• 
Because of the chiral character of the Dirac fermion a, 
this tunneling term raises only forward scattering, the 
effect of which is merely to shift the chemical potential. 
As a result, the chiral Dirac fermion is still gapless. Go- 
ing back to the Majorana fields, we conclude that the 
two chiral Majorana edge modes are stable against suf- 
ficiently dilute impurities. In contrast, in the case of 
the dry-wave pairing, the long-range tunneling via nodal 
excitations exists even in the absence of impurities, as 
depicted in FIGOjF). In this case, an energy gap opens 
around k y ~ 0, and the Majorana mode disappears even 
for a relatively large value of L, for which the Majo- 
rana mode still exists for the d x i_ y i wave pairing, (see 
FIGsfHB) and (E)) 

We, now, consider the stability of the Majorana mode 
in a vortex core. In addition to the localized zero en- 
ergy Majorana solution, there are also delocalized states 
caused by gapless nodal excitations with the finite Fermi 
momenta. When there are multiple vortices in the sys- 
tem under consideration, these delocalized states raise 
long-range tunneling between spatially separated vor- 
tices, which may destroy the zero energy Majorana mode. 
In the system with odd numbers of vortices, one Majo- 
rana mode in a vortex core survives. However, in the 
case with even numbers of vortices, the Majorana mode 
disappears unless they are separated enough from each 
other. 

Topological order and topological protection of Majo- 
rana fermions — The above consideration strongly im- 
plies that there is a topological order which ensures the 
stability of Majorana fermion modes even for nodal su- 
perconductors with bulk gapless excitations. However, 
in sharp contrast to a gapful topological order, the bulk 
Chern number z/ch is not well-defined for our gapless sys- 
tem. Nevertheless, we clarify here that the parity of the 
Chern number (— l)" 011 is well-defined even for nodal su- 



perconductors. The parity of the Chern number ensures 
the stability of the topological order in our system. 

Let us first try to define the Chern number in our gap- 
less system. The simplest way to do this is to introduce 
a small perturbation eliminating all nodes (i.e. gapless 
points) in the spectrum. For instance, adding a small 
id xy term in the gap function, we can easily remove all 
the nodes in our d x i_ y i superconductor. After removing 
the nodal points, the Chern number can be evaluated in 
the standard manner. This procedure, however, does not 
work well after all. The problem is that the value of the 
Chern number depends on the perturbation we choose. 
As a result, one can not have a unique definition of the 
Chern number for gapless systems. 

On the other hand, we find that this procedure does 
define the parity of the Chern number uniquely. From the 
particle-hole symmetry, the parity of the Chern number 
is recast into 
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where Ai(k) is the "gauge field" defined by the bulk band 
wave function \u n (k)), Ai(k) = ^2 n (u n (k)\d ki u n (k)) , 
and Ti is the time-reversal invariant fc-points, I\ = i 2,3,4 — 
(0, 0), (71-, 0), (0, tt), (tt, tt) [31]. Then, for the Hamiltonian 
(Q}, we can show that 



(-1)" 



n 



sg^[4, 



Ar, 



(3) 



i=l,2,3,4 



irrespective of the perturbation (such as id xy term) we 



choose [32l l33j |. This means that we have a unique value 
of the parity in the limit of id xy — > 0; i.e. the parity 
of the Chern number (— l) Ucil is well-defined even for 
nodal superconductors, although the Chern number i^ch 
itself is not. The parity of the Chern number charac- 
terizes the topological phase in nodal superconductors. 
For (— l)^ Ch = — 1, there exists an odd number of topo- 
logically stable Majorana fermion in the edges and in 
a vortex core for nodal superconductors. For example, 
for the model ([T]) with — At — \i < h < — /1, we obtain 
(— l)' ych = —1 from ([3]). Thus, the existence of the gap- 
less Majorana edge mode in FIGQ] is characterized by 
the odd parity (— l) Uch = — 1. On the other hand, for 
(_]_)"ch = 1. there is no topologically stable Majorana 
fermion. We emphasize that the formula (J3j> is appli- 
cable only to systems with particle-hole symmetry, and 
thus, the topological order in gapless systems is specific 
to topological superconducting states. 

In addition to the parity of the Chern number, one 
can consider another topological number dubbed ID Z 2 
invariant [3l|. The ID Z 2 invariant (-1)^°^ is intro- 
duced as a line integral along a specific time-reversal in- 
variant path Cij passing through Ti and Tj . In a similar 
manner above, it is shown that the ID Z 2 invariant is 
well-defined even for our nodal superconductor, and we 



(.1) 



;S gn[e£.+Af,.- 
with 4t — n < 



^sgn^.+A*.-/* 2 ]. 
h < —fi, this formula 



= — 1 for both d x 2_ y 2 and d^j, super- 



obtain (-l)<d c « 
For the model 
yields (-I)"! 01 
conductors. From the bulk-edge correspondence, this Zi 
invariant determines the location of the Majorana edge 
fermions at k y ~ 0, as illustrated in FIGs[T](A),(B), and 
(D). Since the ID Z^ invariant is associated with a local 
structure in the Brillouin zone, its non-triviality does not 
directly lead to the topological stability. However, the ID 
Zi invariant is useful for identifying the location of zero 
energy Majorana edge modes, as shown above. On the 
other hand, the odd parity of the Chern number intro- 
duced above definitely characterizes the global nontrivial 
topology of the Hilbert space, ensuring the topological 
protection mechanism of Majorana modes. 

The above consideration can be straightforwardly gen- 
eralized to a general multi-band nodal superconductor. 
In this case, when the superconducting gap vanishes or 
becomes sufficiently small at the time-reversal invariant 
fe-points r^, the parity of the Chern number is evalu- 
ated as (-l)" Ch = nQ,i=i,2,3,4 s g n [ f a( r i)]i where £ a (k) 
is the normal dispersion of the superconductor. The in- 
dex a specifies an energy band including the spin degrees 
of freedom. Therefore, when the Fermi level is located 
close to odd numbers of time- reversal invariant fe-points, 
and the superconducting gap vanishes at these points be- 
cause of the symmetry requirement, the NCS nodal su- 
perconductor possesses topologically protected Majorana 
fermion modes under an applied small magnetic field. 

Experimental detection of Majorana fermions — For 
the experimental detection of Majorana fermions in 
nodal superconductors, one promising approach is to 
exploit an interferometry measurement proposed for 
a superconductor-topological-insulator junction in refs. 



M 



35| . We consider a setup similar to those proposals, 
but with a difference that, instead of a superconductor- 
topological-insulator junction, a bulk d-wave Rashba su- 
perconductor is used. The contribution from nodal ex- 
citations to the conductance in the d-wave pairing state 
vanishes like ~ T at sufficiently low temperatures, and 
thus, the current is dominated by that carried by two Ma- 
jorana edge modes. The dependence of the conductance 
on the parity of vorticity inside the superconductor sig- 
nifies clearly the Majorana fermion contributions 34Ll35|. 

The non-Abelian nodal superconductor considered 
here is also realizable in an interface between a cen- 
trosymmetric nodal superconductor such as High-T c 
cuprates and a semiconductor, as considered in the 
case of the s-wave pairing state by Sau et al. and 
Alicea la, [la] ■ I n such a system, because of the consider- 
ably large superconducting gap, the experimental detec- 
tion of Majorana modes may be easier. 

Summary — We have demonstrated that even in nodal 
superconductors such as the <i-wave pairing state, Majo- 
rana fermion modes, which are topologically protected 
against weak perturbations and lead to the non-Abelian 



statistics, are realized under a certain realistic condition, 
in spite of the existence of bulk gapless nodal excitations. 
Our results establish a concept of a gapless topological 
phase, and open the possibility of detecting Majorana 
fermions in various NCS superconductors with gap-nodes 
found in real materials. 
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